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Comparisons are made between the av- 
erage magnetic flux density for a three- 
axis circular coil probe and the flux 
density at the center of the probe. The 
results, which arc determined assuming 
a dipolc magnetic field, provide infor- 
mation on the uncertainty associated 
with measurements of magnetic fields 
from some electrical appliances and 
other electrical equipment. The present 
investigation extends an earlier treat- 



ment of the problem, which did not 
consider all orientations of the probe. 
A more comprehensive examination of 
the problem leaves unchanged the con- 
clusions reached previously. 
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1. Introduction 



This paper reconsiders a problem related to the 
measurement of power frequency magnetic fields 
from electrical appliances using three-axis circular 
coil probes. Specifically, it reexamines the differ- 
ences between the average magnetic flux density as 
determined using a magnetic field meter with a 
three-axis circular coil probe and the magnetic flux 
density at the center of the probe, B u , assuming the 
field is produced by a small loop of alternating cur- 
rent, i.e., a magnetic dipole. The "average" arises 
as a consequence of the averaging effects of the 
coil probes over their cross sectional areas when 
placed in a nonuniform magnetic field. The differ- 
ences between the average magnetic field and Bo 
can be regarded as measurement errors because 
the center of the probe is normally considered the 
measurement location. The magnetic dipole field is 
chosen as the relevant field because its geometry 
provides a good approximation of the magnetic 
field produced by many electrical appliances [1]. 



The average magnetic flux density measured by 
a three-axis magnetic field meter, B^, is also re- 
ferred to as the resultant magnetic field and is de- 
fined as [2] 



B mi = VB?+Bf + Bl 



(I) 



where B h B 2 , and flj are average root-mean-square 
(rms) magnetic field components determined by 
each of three orthogonally oriented coil probes. 

Differences between B lv3 and Bo are calculated 
as a function of r/a where r is the distance between 
the magnetic dipole and the center of the probe, 
and a is the radius of the three-axis probe. In addi- 
tion, differences between B„y and B are examined 
for different orientations of the magnetic dipole 
and rotations of the three-axis probe. Because the 
relative orientation of the dipole and three-axis 
probe is not known during most measurement situ- 
ations, there is a distribution of possible differ- 
ences between B av .i and Bo, and these differences 
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collectively represent a source of measurement un- 
certainty for a given rja. What will be of interest in 
this paper is the largest difference that occurs 
between B^ and B as a function of rla (for all 
possible orientations of the dipole). This largest 
difference is designated AB mMi . 

This investigation extends an earlier treatment 
of the problem which considered different orienta- 
tions of the dipole, but not all possible orientations 
of the three-axis probe [3]. The maximum differ- 
ence between B aV 3 and Bo, AB mx o, is found by a nu- 
merical search during which B aV 3 is determined by 
numerical integration. The major advance over the 
earlier study is the development of an expression 
giving the average magnetic flux density for a circu- 
lar coil probe for any position and orientation of 
the probe in the dipole magnetic field. This devel- 
opment allows the search for AB mu .i to consider 
"all" possible rotations of the three-axis probe. 
The extended search is shown to leave unchanged 
the values of ABm^a that were determined by the 
earlier treatment. 



2. Expression for Average Magnetic Field 

In the derivation given below, it is assumed that 
the cross sectional areas of the wire in the coil 
probes and the opposing magnetic fields produced 
by currents induced in the probes are negligible. 
We also assume that the three orthogonally ori- 
ented coils of the three-axis probe have circular 
cross sections of equal area. These assumptions ei- 
ther can be met in practice or can be taken into 
account by a calibration process. 

The average magnetic flux density, B m for a sin- 
gle circular coil probe with cross sectional area^ is 
given by 



B™=\\\ Bn&A, 



(2) 



where cL4 is an element of probe area, n is a unit 
vector perpendicular to A , and B is the magnetic 
flux density. In spherical coordinates, the magnetic 
flux density for a small current loop of radius b is 
[4] 



B = 2*3 cosfor-f . 3 sinffag, 



(3) 



where im, is the permeability of vacuum, / is the 
alternating current, and u r and u e are unit vectors 
in the directions of increasing r and 0, respectively. 
The assumption is made that b<$r, and the sinu- 



soidal time dependence of the field has been sup- 
pressed. The value of Bo is given by the magnitude 
of B [Eq. (3)]. Figure 1 shows the spherical coordi- 
nates r and 0, a small current loop at the origin of 
the coordinate system, and a sketch of a three-axis 
probe. The center of the probe coincides with the 
origin of the prime coordinate system x',y', andz'. 
The coil probes are labelled PI, P2, and P3, have 
unit normal vectors «i, n 2 , and n lt respectively, and 
are shown in Figure 1 (inset) for illustrative pur- 
poses as being in the directions of prime coordi- 
nates. The orientation of the magnetic dipole with 
respect to the position of the probe is characterized 
by the angle 6. 



loop of 

alternating 

current 




* y 



>x° 



Fig. 1. Three-axis magnetic field probe with its center atx=xu, 
y = Q, and 2 =2 . A small current loop producing a dipole mag- 
netic field is located at the origin of the tmprimed coordinate 
system. The unit vectors tti, n 2 , and W3 arc normal to the areas of 
probes PI, P2, and P3, respectively. Changes in the angle B cor- 
respond lo varying the orientation of the dipole with respect to 
the probe. 

For our purposes, it is convenient to express B in 
terms of Cartesian coordinates. The magnetic flux 
density is then [2] 



B=i 



3Q2.;3Q2 + i£J?*L 



2r : 



+J 



2r s 



+ *ia?->). 



(4) 



where r = [x+y+z], i,j, and k are unit vectors 
for the Cartesian coordinates, and C is the constant 
IMJb 2 /2. 
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The goal is to develop an expression for B av at an 
arbitrary point which can be evaluated for any ori- 
entation of the coil probe. The value of B av3 can 
then be found by combining the rms values of S» v 
from three orthogonal directions according to Eq. 
(1). The approach described below for obtaining 
the desired expression for B& is to transform the 
problem into the coordinate system of the coil 
probe. In this coordinate system, the unit vector 
normal to the plane of the coil coincides with the 
"z-axis", B is expressed in terms of the probe coor- 
dinates, and the integration over the area of the 
circular coil probe is carried out numerically in po- 
lar coordinates. 

We begin by considering, without loss of general- 
ity, a three-axis coil probe with its center &tx=x n , 
y = 0, and z=-Za, where x a =-rsin8 and z n =rcos8 
(Fig. 1). We then focus on coil probe PI and its 
unit vector «j after it is rotated through angles a x 
and «2 with respect to the prime coordinate system 
as shown in Fig. 2. The unit vectors n 2 and n 3 will 
also change in orientation to maintain their orthog- 
onal relationship, but are not shown for purposes 
of clarity. In the prime coordinate system, iii is 
given by 

«i=/sinai costa-f/sinaisinai + itcosai. (5) 

By examination, 

n2 = isin(«i + 90°)cosa2-f-jsin(ai -t-90")sina2 
+Jtcos(ff, + 90°) 

- i COSaiCOS«2 +j cos«isina2 — k sinaj. (6) 

The remaining unit vector, n 3 , is given by 



JI3 — Hi X«2 

= -i"sina 2 +,/cosa;2, 



(7) 



but it also can be determined by examination, i.e., 
at is replaced by 90° and ai is replaced by 012 + 90° 
in Eq. (5). For this case, JI3 is constrained to lie in 
the x'-y' plane. Later this constraint is removed. 

The coordinate system of the probe is reached by 
the following transformations: 

(i) translation of the origin to the origin of the 
prime coordinates shown in Fig. 1, 

(ii) rotation of the prime coordinates through 
angle 05 about the 2 '-axis, yielding the double- 
prime coordinatesjc",y", 2", as shown in Fig. 2, and 

(iii) rotation of the double-prime coordinates 
through angle ay about the y"-axis yielding the 
triple-prime coordinates x m ,y m , z m (Fig. 2). In the 
triple-prime coordinate system, the normal vector, 
«j, given by Eq. (5), is along the z "-axis as desired. 



P1 




oy 



Fig. 2. (a) Geometry of unit vector «j and coordinates after 
rotation of the prime coordinates through angle aj about the 
z'-axis and after rotation of double-prime coordinates through 
angle c^ about y"-axis. 

Transformation (i) is given by 



X=Xu+x' 

y=y' 

Z -Za+z'. 



(8) 



The first rotation of coordinates (ii) is given by 

x' =Jc"cosa2-y"sina2 
y ' =x"s\na 2 + v "coses 
z'=z", (9) 

and the second rotation (iii) is given by 

z" ~z '"costri — x "'sinai 
x"=z'"sma l +x"'co$a, 
y"=y'". (10) 

From Eqs. (8)~(10), we have 

x =*b+ (z^'sina) +;t'"cosai)cosa2 —y "'sina 2 
y = (r'"sinai+x'"cosar[)sin«2+y"'cosa 3 
z=zo+2"'cosai-x'"sin«i, (11) 

which, when substituted into Eq. (4), expresses B in 
terms of the probe coordinates. A simplification of 
Eq. (11) is had by noting that the integration over 
the area of the probe [Eq. (2)] occurs in the* m -y ™ 
plane, i.e., z'"=0. Contributions to S av3 from each 
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of the coil probes is found by using the appropriate 
normal vector [Eqs. (5)-(7)] during the integra- 
tions. As noted above, the integration is carried out 
numerically (using a double Simpson's Rule) in po- 
lar coordinates, i.e., 



jc" r = pcos^, 0*Sp=Sa, 0^i!f<2Tr 

y'"=psint/r, 

d4=ck'"dy'" = pdpd^. 



(12) 



The accuracy of the numerical integrations was 
checked by increasing the number of divisions be- 
tween the limits of integration for p and $. The 
results reported below were not affected by further 
refinements of the intervals used during the inte- 
grations. 

In the search for AB m „3, it will be necessary to 
perform rotations of the three-axis probe about the 
z m -axis or m direction (see Search Protocol below), 
i.e., the unit vectors ni and rts for probes P2 and P3 
are rotated about m. This removes the constraint 
noted earlier that n 3 lies in the x'-y' plane. Because 
the integrand for B m is in terms of the a angles, 
relationships must be found between the angle of 
rotation about the z"'-axis, designated as $, and 
the a values that appear in the integrands for P2 
and P3. These relationships are found by examin- 
ing the unit vectors for the probes n 2 and /i 3 as they 
rotate about the z "-axis or «i direction. 

In Figs. 3 and 4, consider n t in a direction char- 
acterized by the angles a w and a w , i.e., 

nt=isinorioCOSOM+jsinci;iosinor2fl + ^cosain. (13) 

From Eqs. (6) and (7), we have 

tit = (sin(aio + 90°)cosa2o +jsi'n(aio+ 90°)sinaa) 
+*rcos(ff 10 +90°), (14) 

and 

/la^/sinf^costa^-l- 90°) +jsin(90 ,) )sin(a 31 + 90 o ) 
+ fccos(90°). (15) 

The trigonometric expressions in Eqs. (14) and 
(15) are not simplified in order to aid the reader in 
seeing the relationships between the three unit vec- 
tors. 

Following a counterclockwise rotation of 4> de- 
grees about the z^-axis, the ct's will increase in 
value and the et\% will decrease in value in the ex- 
pressions for iti and « 3 . These changes also occur in 
the expression for the magnetic flux density B. 




n a <<*-0) 



n 2 (<» * 0) 



"?(^ = 0) 



Fig. 3. Geometry of coordinates and unit vectors after Unit vec- 
tor n 2 is rotated <j> degrees about n, or z '"-axis. The rotation of 
Hi is not shown for purposes of clarity (sec Fig. 4), 




" 3 (* = 0) 



circle of 

rotation 
(radius c) 



'"*(^0) 



"j<S> = 0) 



Fig. 4. Geometry of coordinates and unit vectors after unit vec- 
tors n 3 and « 2 are rotated through angle <j> about n t or z "-axis. 
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After a rotation of 4> degrees, a line along the 
unit vector ni will intersect the circle of rotation in 
the x m -y" plane at a point given by (Fig. 3) 



x m =ccos<{> 
y*"=csin<£ 



(16) 



where the radius for the rotation has been arbitrar- 
ily taken to be some constant c. From Eqs. (10) 
and (16), the same point in the double prime coor- 
dinate system is 



z"~ -jf"'sincrio= -ccos^>sinaio 
x" —x "cosaio = ccos^coscrio 
v"=y™=csin^>. 



(17) 



The increment to q*b for unit vector rti after rota- 
tion <£, Ba (Fig. 3), can be found from the expres- 
sion for its tangent, i.e., 

tan&= 4 =J an£ 
X cosaio 

fe-tan-(-^Y (18) 

\ COSOki/ v ' 

Prior to the rotation, the angle with respect to the 
^'-axis for the unit vector m is given by oio + 90 D 
(Fig, 3). After the rotation, the corresponding an- 
gle will be 5i2+90° where 5] 2 <a I0 . The value of 8a 
is found by noting that after the rotation <f>, the line 
/ from the origin to the projection of the circle of 
rotation onto the x"-y* plane is given by 



/=V(*") 2 +(y") 2 , 



(19) 



and that the tangent for $ l2 is just [z7/| (Fig-3). 
From Eqs. (17) and (19), 



Iz 



= tan* 



cos^sinoiu 



* VcosVcos 2 aio + sin 



5) (20 » 



Thus, following a rotation of 4> degrees about n h 
aw and aio + 90° will be replaced by 020 + fe and 
5i2+90 o , respectively, in the expressions for »2 and 
B during the calculation of B K for probe P2. 

The or values for m can be determined with a 
similar analysis. Following a rotation of $ degrees 
about the z™-axis {Fig. 4), a line along the unit 
vector «3 will intersect the circle of rotation in the 
x m -y" plane at the point 



x** = ~csin<£ 
y m = ccos<f> 
z"=0. 



(21) 



From Eqs. (10) and (21), the same point in the 
double prime coordinate system is 



z" = csin<£sinaiti 
x H — — csin^cosaio 
y"=ccas<f>. 



(22) 



Following rotation <j>, the angle 020 + 90° for nj 
[Eq. (15)] will increase by an amount S21 as shown 
in Fig. 4. The increment, fe, can be found from the 
expression for the absolute value of its tangent, i.e., 



tanS B = 



= tan 4> cosaio, 



523 = tan" 1 (tan^cosaio). 



(23) 



Prior to the rotation, the unit vector n% makes an 
angle of 90* with respect to the z'-axis (Fig. 4). 
After the rotation, this angle will decrease by an 
amount 5i3. The value of 5n is found by noting that 
after the rotation <f>, the line m from the origin to 
the projection of the circle of rotation onto the 
x"-y" plane is given by 



m =V(xy + (y"f, 



(24) 



and that the tangent for S n is just z"/m (Fig. 4). 
From Eqs. (22) and (24) we have 



tan8i3=— = 
5u = tan"' 



sin^sinom 



m Vsin 2 0cos 2 aio+cos 2 <£ 

I sin<ftsina-in \ 

Wsin 2 0cos 2 ario+ COS 2 <£>' ' 



(25) 



Thus, following a rotation of $ degrees about «i, 
020 + 90° and 90° are replaced by O20 + 90 o + 523 and 
90°-fiu, respectively in the expressions for m [Eq. 
(15)] and B during calculation of S av for coil probe 
P3. 



3. Search Protocol 

The search for the largest difference between 
5 av j and B 0t ABmixi, for a given distance r from the 
dipole proceeds as follows: 

(i) For a fixed distance r away from the dipole, 
and with & = 01 = 05 = 0, the three-axis probe is ro- 
tated about the z^-axis or iti direction in 2° steps 
(i.e., 4> is incremented in 2° steps). Then B m for 
each coil probe is evaluated and combined accord- 
ing to Eq. (1) for each value of <f> to obtain B, v j, B m3 
is compared with #0, and the largest difference is 
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saved. Because of the symmetry of the problem, a 
total rotation of 90° is required to cover all the 
cases (with 2° increments). 

(ii) The angle «i is advanced in 5° steps and the 
above comparisons are repeated as the probe is ro- 
tated about the z m -axis or m direction. The maxi- 
mum value of «], without duplication of results is 
90°. 

(iii) For each value of a ir ai is incremented from 
0° in steps of 5° and the above comparisons are 
repeated. Because of symmetry arguments, a total 
rotation of 180° is required to consider all the cases 
without duplication. 

(iv) Following the above calculations, different 
orientations of the magnetic dipole are considered 
by changing the angle 6 in 15° steps and repeating 
steps (i) through (iii). The choices of increments 
indicated above were found to provide adequate 
sensitivity for determining AS„ 3 . 

(v) Steps (i) through (iv) are repeated for differ- 
ent values of r. 

A diagram schematically indicating several posi- 
tions for m, and rotations about m, as the above 
protocol was carried out for a fixed value of r is 
shown in Fig. 5. 




Fig. 5. Trajectories of «i during search for &Bw X . The unil vec- 
tors tt 2 and Hi are not shown but maintain their orthogonal rela- 
tionship wilh n throughout the search. 



4. Results and Discussion 

As already noted, an earlier search for ABn^ [3] 
was not as comprehensive as the one described in 
this paper. While the ratio r/a and could be varied 
without restriction, the rotations of the three-axis 
probe were limited to simple rotations about the 
* '-■> y '-> or z '-axis. That is, it was not possible to con- 
sider differences between B 3Vi and B for combina- 
tions of rotations about two or three axes. This 
problem has been overcome with a more general- 
ized expression for B„ v compared to the ones used in 
the earlier calculation. What is perhaps surprising, 
however, is that the AB msi3 values obtained with the 
more comprehensive search protocol are the same 
as previously calculated. That is, ABm^a is negative 
and occurs for all r/a values when 6 = 90°, as previ- 
ously found, and correspond to the AJ? mM 3 values 
determined earlier following simple rotations about 
they '-axis (referred to as "a rotations" in Ref. [3]). 
Numerical values of Afl raajl 3 are provided in Table 1 
as a function of r/a . 

Table 1. Values of AS„„,j as a function of normalized distance 
r/a from magnetic dipole 



r/a 



AS „„,,(%} 



3 
4 
5 
6 

7 

8 
9 
10 
11 
12 
13 
14 
15 



-19.6 
-10.8 
-6.9 
4,8 
-3.5 
-2.7 
-2.1 
-1.7 
-1.4 
-1.2 
-1.0 
-0.9 
-0.8 



5. Conclusions 

The present calculations have determined the 
largest differences between the resultant magnetic 
field, B aV 3, and the field value at the center of the 
probe Ba, assuming a dipole magnetic field. These 
largest differences, designated Afl mM 3, are reported 
in Table 1 as a function of normalized distance, 
r/a, from the center of the dipole and agree with 
values previously found after a far less comprehen- 
sive search [3]. The quantity, AB max j, can be re- 
garded as the largest error due to instrumental 
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averaging effects. As noted earlier, because the rel- 
ative orientations of the dipole and three-axis 
probe are not known for a given r/a under typical 
measurement conditions, there will be a range of 
possible differences between B„ 3 and £J , Thus, ide- 
ally, it would be desirable to determine the distri- 
bution of differences between B^i and fl and treat 
the problem using a statistical approach, but that 
has been left to a future calculation. 

Because the dipole field is a good approximation 
of fields produced by many electrical appliances, 
the information in Table 1 should be taken into 
account when total uncertainties are being deter- 
mined during measurements of magnetic fields 
from appliances. For example, if the resultant mag- 
netic field is to be measured at a distance r from an 
appliance with a combined relative standard uncer- 
tainty [5] of less than ± 10%, magnetic field meters 
with three-axis probes having radii a such that 
r/a ~3 should be considered unsuitable. Three-axis 
probes having radii such that r/a =5 would conser- 
vatively be considered suitable if the combined rel- 
ative standard uncertainty from all other sources 
(e.g., calibration process, frequency response) 
amounted to about 3% or less, since 6.9% -t- 
3.0% = 9.9%, where 6.9% is taken from Table 1 for 
r/a =5. 
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